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Measuring thermodynamic quantities can be easy or not, depending on the system that is being studied. For 
a macroscopic object, measuring temperatures can be as simple as measuring how much a column of mercury 
rises when in contact with the object. At the small scale of quantum electromechanical systems, such simple 
methods are not available and invariably detection processes disturb the system state. Here we propose a method 
for measuring the temperature on a suspended semiconductor membrane clamped at both ends. In this method, 
the membrane is mediating a capacitive coupling between two transmission line resonators (TLR). The first 
TLR has a strong dispersion, that is, its decaying rate is larger than its drive, and its role is to pump in a pulsed 
way the interaction between the membrane and the second TLR. By averaging the pulsed measurements of the 
quadrature of the second TLR we show how the temperature of the membrane can be determined. Moreover the 
statistical description of the state of the membrane, which is directly accessed in this approach is significantly 
improved by the addition of a Josephson Junction coupled to the second TLR. 
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I. INTRODUCTION 


Electromechanical systems are devices which couple mechanical displacement and electrostatic interactions. Measuring phys¬ 
ical properties of such a device at macroscopic scales is relatively easy - Coulomb in his famous torsion balance attached charged 
metal spheres to rods and threads and visually measured the torsion produced by the electrostatic interaction between the spheres. 
The measurement of the torsion allowed him to determine the force acting on the spheres IHGl. At the nanoscale however Q, 
the movement of a nanoelectromechanical system (NEMS) cannot be observed directly. A further complication is that at those 
scales a quantum description of the system is invariably necessary. 

The specific NEMS we are interested in is a suspended semiconductor membrane clamped at both ends, which will be oscil¬ 
lating due to coupling to the thermal modes of the clamps. The smaller is the mechanical element of the NEMS, the stronger is 
the coupling to the thermal modes of the reservoirs that clamp them at the extremities, and usually some enginnered structures 
must be employed in order to decrease its effects 0|. This oscillator can be coupled electrostatically to other devices, allowing 
the transduction of the movement as electric signals (See Ref. (5) for a general review). Previously schemes to measure the 
quadrature phase amplitude |6] and to observe the quantum of thermal conductance (7) have been proposed. It is particularly 
relevant that the detection of the NEMS movement can give direct access to its temperature, a fundamental physical quantity 
GOES' For the NEMS temperature measurement, usually the area under the noise power spectrum of the displacement amplitude 
transduced signal is used as it gives directly the mean phonon number in the steady state (See SHE) for example). However, 
one could ask on how to access the temperature of the mechanical resonator by means of a non-demolition detection scheme. 
Indeed, some previous discussion on non-demolition detection in the accessment of temperature has been given ifTTIfTTll . Also 
a discussion on back-action-evading detection and squeezing generation in mechanical resonators through coupling to radiation 
field were given in fl3l as well as the outstanding implementation for movement detection given in fl4l . 

In this work we describe a circuit in which we can perform repeated quantum nondemolition measurement sof the expected 
value of the number of phonons of the NEMS, and show an analysis of the statistics of the number of phonons in this circuit. 
Through this analysis we show how the temperature of the semiconductor membrane can be directly accessed. Finally, we show 
that by discretizing the measurement of the expected value we can increase its accuracy in a pulsed linear radiation detection 
process Dana. 



FIG. 1. Diagram for a pulsed quantum nondemolition measurement. The drive “D” pulses the interaction between the system “S” and the 
measuring device “M”. The noise reducer “NR” lowers the fluctuations of the meter in each measurement. This cycle, is repeated many times, 
always following the order: state preparation —>• coherent control —quantum measurement. 


In order to perform a measurement on a small oscillator, in a quantum regime |l] [2, 14, 171], a quantum nondemolition 
measurement is necessary. In this type of measurement, the coupling with the meter is such that it does not disturbs the non¬ 
demolition observable of interest. Also, because these measurements are probabilistic by nature, it is also important to improve 
the precision of the meter so that at the end of many measurements there is a significant accuracy in the calculated average. This 
cyclic sequence of non continuous measurements in time is shown on the diagram of Fig. [I] 
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II. RESULTS 

We consider the coupling of two Transmission Line Resonators (TLRs) mediated by a mechanical resonator as considered in 
Ref. (HI and represented in Fig. [2] By assuming the regime of rapid mechanical oscillations at the GHz scale, and that the 
TLR-1, being in a undepleted regime, is treated classically, the interaction Hamiltonian between the NEMS and the field in the 
TLR-2 is given by 


Hi = — ha(t)b^b(a + a^), ( 1 ) 

where a(t) is proportional to the microwave amplitude in the TLR-1 [ 18] and a(a^) is the annihilation(creation) operator of the 
normal mode in the TLR-2. b(b t) is the mechanical resonator annihilation(creation) operator. As one can see in this interaction 
Hamiltonian the phonon number operator b^b is a nondemolition variable. In addition the field in this remaining TLR is affected 
by the presence of a Josephson Junction (See Fig. [2]), whose effect (as demonstrated in the Appendix) is to induce quadrature 
squeezing lfl9ti25l through a parametric term 


H p = -ih^(a 2 -a t 2 ). (2) 

Without loss of generality we take 7 = /3x ^/2 as being real and positive, where is the nonlinear susceptibility and /3 is 
the amplitude the coherent pulse source in the TLR-2. We shall consider a pulsed measurement scenario lfl5l 16] where the 
interaction in Hj is turned on and off rapidly. Note that the unconditional phonon number operator is a quantum non-demolition 
variable. Moreover given that [tfb, H p + Hi] = 0, back-action evading is guaranteed and therefore tfb can be repeteadly 
measured without the back-action noise [26]. Therefore the phonon number statistics is not changed from pulse to pulse. The 
detailed process of this pulsed measurement is described in what follows. For a detailed description of the elements involved 
forlinear detection please see Refs. IT5lfT61. 

The experiment is carried out according to Fig. [2] The Josephson Junction acts as a source of coherent pulses, with time 
interval £ 1 , on the TLR-2 In sequence, a coherent pulse is applied to TLR-1, with amplitude aq and time interval £ 2 . The 
radiation field there is treated classically, as a result of the undepleted regime COD. After these pulses, the density operator is 
given by 

p(r = fi + h) = p (0) e iH ^e iHlt2 . (3) 

Initially, the TLR-2 is prepared in a vacuum state, while the NEMS, which oscillates due to thermal excitation, in equilibrium is 
in a thermal state with mean phonon number N, 


p(0) = ^P(n)|n)(n| 6 ®|0)(0| o , (4) 

n =0 

where P(n) = N n /(N + l ) n+1 is the thermal phonon number distribution in the NEMS, being N = (exp (hv/kBT) — l ) -1 
its thermal number. Given the very shorts pulses with time intervals r = t\ + £ 2 , the evolution can be approximated as above to 
give 


00 

P(t) = - P ( n ) I”) Hi, ® l«n, 7) («n, 7l a > ( 5 ) 

n=0 

with a n (r) = in f * 2 a(t)dt = inA , where A is the area of the interaction pulse, and |a n ,7 ) a is a coherent squeezed state 
given by \a n ,x) a = V(a n )S( 7) \n, 0), where T>(a n ) = exp [inA(a + a^)] is the displacement operator conditioned on the 
mechanical resonator phonon excitation number and 5(7) = exp [(7£i/2)(a^ 2 — a 2 )] is the squeezing operator. 

At the end of a pulse the reduced cavity field state is described by a mixture of squeezed coherent states, which are distributed 
according to the mechanical resonator thermal phonon number distribution P(n). The field in TLR-2 is dumped out through 
the dual IQ mixer measurement scheme of Ref. insmn, detailed in Fig. 0 - In this scheme, the TLR-2 output is directed 
to a microwave beam splitter (a hybrid coupler). The two outputs are then amplified and run through separate IQ mixers. 
The four output currents from the IQ mixers can be correlated in various ways after proceeding to linear detectors fI5l . One 
particular cross correlation gives access directly to all the normally ordered moments of the TLR-2 field at the initial time of 
the measurement. At the end of each pulse a measurement is performed on the output of the second TLR through complexes 
envelopes S c = —a(r) and Sd = ia(r). It is important to remark that the criterion for back-action-evasion for b^b, mainly 
[bU), H p + H r ] = 0, is obeyed in every detection step. 

These pulsed measurements allow that the cross-correlations be computed after all the field has been detected and the TLR-2 
returned to the vacuum state ready for the next pulse. For each individual pulsed measurement the NEMS is to be found on an 
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(d) 


FIG. 2. In the proposed experiment, the capacitive coupling between two TLRs is mediated by the vibration of a NEMS. The TLR-2, on the 
right is directly coupled with a Josephson Junction (JJ) for noise reduction. The pulsed measurement sequence is given as follows: (a) Step 
one. Initial state preparation p(0); (b) Step two. Application of a pulse in TLR-2 by a coherent pulse source of generating p(ti); (c) Step three. 
Application of a pulse in TLR-1 by another coherent pulse source generating p(r); (d) Step four. The output of the TLR-2 is directed to a 
microwave beam splitter (a hybrid). The two outputs are then amplified and run through separate IQ mixers. The four output currents from 
two IQ mixers can be correlated in various ways. The particlar correlation employed here gives access directly to all the normally ordered 
moments of the cavity field at the start of the measurement. The cross correlations are computed after all the field has been detected and the 
TLR-2 returned to the vacuum state ready for the next pulse. (Note: the pulse applied in step 3 does not interact with the JJ, for not producing 
sufficiently large charging energy.) 
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specific phonon number excitation m according to the thermal distribution P(m). The post-selected state of the TLR-2 field 
conditioned on this instantaneous (at the time of measurement r) NEMS excitation m is given by 


Tr b {n m p(r)n TO } 
Tr ab {n m p(r)n m }’ 


( 6 ) 


where H m = \m) (m\ b , p(r) is given by Eq. <[5j», and Tr ab {II m/ o(r)II m } = P{m) is the probability to obtain m. Therefore the 
expected value of the phase quadrature in each measurement 


00 (m) ( T ) = \ot m ,i) , 

has a variance 


(7) 


(( Ay ) 2 > (m) W= e “ 27tl - (8) 

The TLR-2 field pre-selected state p a (r) = P(m)pa (t) average over all the quadrature measurements as (Y) (r) = 
P( m ) {Y)(m) ( r ) P rov iding the phonon number (and consequently the temperature) measurement. 

The computation of the total variance on F depends on the variance of the thermal NEMS by adding the variance of each 
individual measurement, and is dependent on the distinguishability of each individual detection. The squeezing enables the dis¬ 
junction of the individual conditioned TLR-2 field detections, allowing to attribute the correct excitation m influencing p^ (r), 
as depicted in the phase space in Fig. [3] Therefore the squeezing improves the statistical resolution, being necessary that 
jti > —ln[y/2A] for that to be ensured. The ellipses are centered at n2A , n = 0 , 1 , 2 ,... , and their minor and major axes 
represent the fluctuation of the measurement of each quadrature. Without the squeezing an incorrect attribution could be given 
to m through the fluctuation ((AF) 2 )^ (r) = e _27tl . 



FIG. 3. Phase space of the TLR-2 in t — r for 7 t\ > —ln[y/2A]. The minor (major) axis of the ellipses represent ((AF) 2 )^ (r) =§ e 27tl 

(((AX) 2 )^ (r) = e 27tl ), the variance of each measurement of the quadrature F ( X ). The separation of the ellipses along F is equal to 
pulse the interaction area 2 A. 


The resulting average over many pulsed measurements is given by the quadratures 

(Y) (r) = Tr {* (a 1 - a) p(r)} = 2AN, (9) 

and 


(X) (r) = Tr{(a + a f ) p(r)} = 0. (10) 

Therefore, by averaging over many identical pulses we can reconstruct the mean phonon number of the NEMS, N , and thus 
deduce its temperature, T = hv/[k,B In (A -1 + 1)] . In a similar fashion we could also measure two-time correlation functions 
for the mechanical resonator. Instead here we focus on the quadratures variance and the increasing of accuracy in determining 
the phonon statistics of the mechanical resonator due to the presence of the squeezing term in Eq. ([2]). 

Through the pulsed measurement the variance of the quadrature F is given by 


<(Ay) 2 ) = 4A 2 N (TV + 1) + <(Ay) 2 > (m) (r), 


( 11 ) 


being ((AF) 2 )^ (r), the variance of each individual measurement given by Eq. | 8 j). N(N + 1) appearing in Eq. (11) is the 
variance of the thermal distribution of the NEMS. The squeezing induced by the Josephson Junction allows a reduction of the 
noise on each measurement, in contrast with the situation without the Josephson Junction, for which ((AF) 2 )^ (r) = 1. In 

the limit of large squeezing the relative uncertainty AF = 7 /((AF) 2 )/ (F) is approximately given by 


AF « y/l + N~\ 


( 12 ) 
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being independent on the pulse area and reapidly reaching the threshold 1 with increasing N. 

Extending the previous procedure for measurement of the TLR-2 field ordered moments we can reconstruct the Wigner 
quasiprobability distribution [27, 281] 


m«) = £ / 


d z A 


<(qt)"0(-A*rA\ _i |A | : 


+aA* —ck * A 


(13) 


for the field state, giving 


W(a) = £ 

n 


p — ^Re 2 (a)e~ 2 ' ytl — ^ [Im(a)—a n ] 2 e 2ytl 

27T 


(14) 


The signature of the NEMS phonon distribution P(n) is clearly imprinted in the Wigner distribution profile in Fig.Q and can 
accessed through the marginal distibution V = / dRe(a)W(a). 



(a) 



(b) 

FIG. 4. (a) Wignerdistibution for the TLR-2 field, as given in Eq. for N = 1, A = 1 and e 2ltl = 50. Integrating the Wigner distribution 
with respect to variable Re(a), V = f dRe(a)W(a ), we obtain the visualization of the NEMS thermal phonon distribuiton P(n) depicted 
1291 in the histogram (b). 
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III. DISCUSSION 

Detection of thermodynamical quantities associated to nanoscopic movement is important in many aspects. Particularly, 
the detection of the corresponding temperature of such tiny devices is relevant for physical characterization Eiia and latter 
usage in further applications. Here we have proposed a mechanism to measure the average number of phonons of a NEMS 
through a pulsed non-demolition detection scheme. Since here the NEMS phonons are due to thermal excitation, a direct 
way to characterize the temperature of the device as well as it statistical properties in a non-demolition way is derived. The 
proposed scheme to measure the temperature of the NEMS is experimentally feasible with nowadays technology. The interaction 
between semiconducting NEMS 0] 111 El El ES [311 and superconducting TLR |'32~34] has already been carried out 151 351-381. 
Experiments with The IQ mixer measuring method were also performed fT 6 l . Therefore we expect that the presented scheme be 
readily implemented. The mechanical resonator is assumed to be at thermal equilibrium at the measurement stage, and therefore 
dissipative effects, which are strong for those devices play a key role for reaching this equilibrium. That is to be reached prior 
the sequence of pulses of the detection scheme, otherwise it would be expected a continuous decrease of the phonon number, 
which would be typical in a transient regime. The more relevant dissipative effect comes in fact from the superconducting charge 
qubit decayng time l39l (See m for a compreensive review of some recent experimental numbers). At the frequencies of GHz 
necessary for the present porposal the charge qubit decay will affect the squeezing of the TLR 2 field, and therefore will affect 
the accuracy of the detection process. However given the ability to perform short and fast pulses in those devices IT31 we do not 
expect it to be detrimental for the present proposal. A more detailed investigation on that is required. 
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Appendix: Derivation of Hamiltonian {2}. 

We consider the model in figure ([5]), where the three lowest Josephson Junction charge states form a three-level artificial atom 
is in a configuration where the ground (| g)) and excited (|e)) states are separated by an intermediary atomic level (|z)). The 
modes with frequency Qi, D 2 and D 3 are coupled to the transitions | g) ^ |i), | i) \e) and | g) |e) with coupling constants 
gi, g 2 and g 3 , respectively. 

The Hamiltonian in the rotating wave approximation is 

H = H 0 + V, (A.l) 

being 

3 

Hq = hVtj a j a j + m e a ee + + hfl g a gg , (A.2) 

3 = 1 

and 

v = h ( gia\(Tig + g 2 a\a e i + g 3 a\a eg + h.c.'j , (A.3) 

where a^ = | j) (k\ is the atomic transition for j, k = g,i, e. Writing H in the interaction picture, by applying Uo(t) = e ~ lHot / h , 
and transforming it through U\ = e -^t(a ee +a gg ) ? we 0 bt a i n 

— = -A (a ee + dgg) 

+gi a lVgi + 92^2^ie + g3 a 3°’ge 

+9l a \&ig + 92 a 2 a ei + 93 a 3 a eg- (A.4) 

Therefore, the Heisenberg equations of motion for the coherences cr g i and <r ie are 

+ gidi (on ®gg) 3 “ 93® j 3&ie 92 ® J 2 ( ^eg 5 


(A.5) 
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(a) 



l*‘> 


FIG. 5. (a) Scheme for generation of a parametric oscil l ator, v ia the pulse applied to the TLR-2, coupled to an artificial three-level atom; (b) 
Energy-level diagram, observing the Hamiltonian ( |A.2| ), \A.3\ and A = —(Qi — Q g — Qi) = Q e — Qi — Q 2 (with Q g = —uj, 9l e = uj and 
Qi — —cji) and Q 3 = Vl e — Vt g — S. 


and 


ter j. 


-Acr ifi 


920-2 


He 


s) + S3 a : 


3®ig 91&1& ge- 


(A. 6 ) 


Assuming that initially the artificail atom is prepared in the state | i) and that both A and S and large enough so that the states \e) 
and \g) are not significantly populated so that an = 1 , in the limit where A >> |#i | , | # 2 1, | # 3 1 , S, we have the adiabatic solution 
with dig = &i e = 0. Next, the expression of a g i and an, we can consider A 2 — g 2 a\a 3 = A 2 , the terms with a gg , cr ee , a ge and 
a eg can be neglected, resulting in the solution 


^ 9 1 . #2#3 f 

^ + ^2“ a 2 a 3 


(A.7) 


and 


^ - 92 + , 9i 93 n + 

a ie = + ^“ a l a 3* 


(A. 8 ) 


In this way, preparing the artificial atom in intermediate state | i), and replacing (A.7) and (A. 8 ) expressions in Hamiltonian 
(|A.4|), we obtain the effective Hamiltonian of the form 


Heff _ 9i 

2 h A 


a\ai + 


92 t 

g4a 2 + 


919293 

A 2 


aia 2 a\ + 


9 i 929 s f t 
A 2 a|a^a 3 , 


(A.9) 


being g\ and g 2 real. Now, applying the transformation U 2 = e ^(sM a r+# 2 a 2 a 2 )/ A ? a3 — being a coherent pulse source of 
amplitude [3 , and #3 -A Q 3 e~ l5t so that its frequency is adjusted to 5 = (gf + g 2 )/A, we obtain 

H p = -ih^(a 2 -a t2 ), (A.10) 

with ai = a 2 = a, = S7 2 = cc, 7 = 4gig 2 (/ 3 /3/A 2 . Note that in this case, the effective frequency is satisfied for 

n 3 = 2u - {gl + g 2 )/ A. 
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